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Table 1: Update of the ensembles produced with j3 = 3.9 by the ETM collaboration. For the other j3-values 
see table 1 of Ref. [||]. We give the lattice volume I? x T, the twisted mass parameter a^lq, the hopping 
parameter (C = 1 /(8 + lam^) and the trajectory length T. In addition we provide values for the integrated 
autocorrelation time of two typical quantities, the plaquette P and the pseudo scalar mass amps, in units of 
T = 0.5. 



1. Introduction 

We present an update of the «/ = 2 results obtained by the European Twisted Mass collabora- 
tion (ETMC) for the pseudo scalar mass and decay constant. A good understanding of continuum, 
thermodynamic and chiral limits is essential in order to obtain reliable results, which can eventually 
be compared to experiment. The physics of the Ught pseudo scalar meson is a prime example for 
investigating these extrapolations, because its mass and decay constant can be obtained with high 
precision in lattice simulations and chiral perturbation theory (;tPT) is best understood for those 
two quantities. In addition, such an investigation allows to extract other quantities of phenomeno- 
logical interest, such as low energy constants and quark masses. 

First results for the pseudo scalar mass and decay constant obtained for the large scale simu- 
lations of the ETM collaboration can be found in Refs. lll]-@]. 

ETMC has generated large sets of gauge configurations for different values of the coupling 
constant (jS = 3.8, a ~ 0.1 fm; jS = 3.9, a ~ 0.085 fm; jS = 4.05, a ~ 0.065 fm), for various 
volumes (2.1—2.8 fm) and a number of bare quark masses corresponding to pseudo scalar meson 
masses ranging from ~ 260 to ~ 700 MeV. The list of ensembles at j3 = 3.9 can be found in 
table |l], which contains the newly generated ensemble S7 corresponding to a pseudo scalar meson 
mass of about mps ~ 265 MeV. For the other j8 -values we refer to table 1 of Ref. [Q]. 

In the gauge sector we employ the so-called tree-level Symanzik improved gauge action (tl- 
Sym) [^. The fermionic action for two flavours of maximally twisted, mass degenerate quarks in 
the so called twisted basis [^, 0] reads 

Stm = [D[U]+mo + wA x{x)] , (1.1) 

X 

where mo is the untwisted bare quark mass tuned to its critical value mait, piq is the bare twisted 
quark mass, is the third Pauli matrix acting in flavour space and D\U] is the Wilson-Dirac 
operator. 

At maximal twist, i.e. mo = nicnu physical observables are automatically ff{a) improved with- 
out the need to determine any action or operator specific improvement coefficients (for a re- 
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Table 2: Comparison of some selected quantities for which an isospin splitting can occur for twisted mass 
fermions. Rq denotes the measured relative size of the splitting. 



view see Ref. With this being the main advantage, one drawback of maximally twisted mass 
fermions is that flavour symmetry is broken explicitly at finite value of the lattice spacing, which 
amounts to 0'{a^) effects in physical observables, as will be discussed later Note that in the fol- 
lowing we shall refer to the charged pseudo scalar meson mass as mps or m^^ and to the neutral 



one as Wpg. 



For details on the set-up, tuning to maximal twist and the analysis methods of the ETM col- 
laboration we refer to Refs. ||^, ^,|^. Recent results for light quark masses and decay constants, the 
light baryon spectrum and the r]' meson aie available in Refs. and Ref. [11], respectively. 

We shall only consider the ensembles at j3 = 3.90 (B-ensembles) and j3 = 4.05 (C-ensembles, see 
table 1 of Ref. |@]) in this proceeding contribution, because tuning to maximal twist at j8 = 3.8 was 
not sufficiently accurate at the lowest quark mass values for the observables considered here. 



2. Results 



Flavour Breaking Effects 



Flavour breaking effects have been investigated by ETMC for several quantities. In figure |l(a) 
we plot ?'o(('Mp5)^ — (/Mps)^) as a function of {a/ro)^. It is visible that mass splitting of the charged 
to neutral pseudo scalar meson is large. However, the measured splittings are compatible with 
being an ^(a^) effect, as expected, and they vanish towards the continuum limit. 

All other possible splittings investigated so far are compatible with zero. In table § we have 
compiled the relative difference Rq = {0 — 0')/0 for some selected simulation points and ob- 
servables O. Here O (O') denotes the charged (neutral) quantity for mesons and A+ (A++) for 
baryons. The values of Rq are well compatible with zero for all observables O besides the pion 
mass. However, some quantities, like the vector meson decay constant /y, are rather noisy, making 
definite conclusions difficult. These results are compatible with a theoretical investigation using 



the Symanzik effective Langrangian 1 12]. 



SU(2) xPT Fits 

As the details of finite size corrections for nips and fys were discussed in Ref. [0] and the issue 
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of continuum extrapolation in fixed and finite volume for selected quantities in Refs. 
shall summarise here only the main results: 



i, we 



1. finite size effects in /ps and mps can be described using chiral perturbation theory in the form 



of the resummed Liischer formula as described in Ref. Q13|]. We denote the corresponding 
correction factors with K'j:^, which depend among others on the low energy constants Ai_4. 

2. within our current statistical precision lattice artifacts appear to be negligible, in particular 
for w^ps,/ps^ the quantities we consider here. In fixed volume and at fixed value of rQtnps 
the results for ro/ps at j8 = 3.9 and j8 = 4.05 are always compatible with each other within 
our small statistical errors, and hence a constant extrapolation to the continuum limit seems 
justified. When we include also a linear term in {a/roY in the extrapolation, a least square 
fit determines a value for the slope that is zero within errors. However, we include this effect 
into our systematic uncertainties. 

We shall now present the results for a combined chiral, thermodynamic and continuum extrapola- 
tion of wips and /ps for the two j3-values j8 = 3.9 and j8 = 4.05. What we present here will extend 
the results given in Refs. ||] by incorporating chirally extrapolated data for the renormalisation 
constant Zp and the Sommer parameter vq/q into the fit. Details on the computation of Zp (using 
the RIMOM) and ro /a can be found in Refs. @, g]. 

We perform combined fits to our data for /ps, mps, r^/a and Zp at the two values of j3 with the 
formulae: 



'"o/ps = '"o/o 



l-2^1og^^j+D,,,(a/ro)^ + r/^L° 



(romps)' = Z^ir^ [1 + ^ log ( ^ j + A„ps {a/roY + 



nNNLO 



(2.1) 
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Figure 2: (a) Data for {ram^^Y / kq^r as a function of r^jXr and Fits A and C. (b) Data for ro/ps as a function 
of r^pLR and resulting curves of Fit B. The vertical lines indicate the fit range. 



with = 2Bol^R/{4Kfo)^ = IBqIXr ,ixr = IX^/Z^Jq = ^/IFq 

NNLO terms [|^, which depend on Ai_4 and Rm and kf, and K'^j^{L) the finite size correc 



jNNLO denote the continuum 



tions [13]. Based on the form of the Symanzik expansion in the small quark mass region, we 



parametrise in eq. (2.1) the leading cut-off effects by the two coefficients Dfy^^^ps- 

At NLO, i.e. setting T^^^'^ = 0, and neglecting finite size corrections for the r 
are the following free parameters to be fitted to the data for a/ps, amps, ro/a and Zp: 

ro/o, roBo, roA^, roA4, {ro/a}p, {Zp}p, D^p^, D^^, 



where we indicate with the notation {■■.}p that there is one parameter for each value a. 

Finite size effects are corrected for by using the asymptotic formulae from CDH, which is 
consistently included in the fit. However, depend on more parameters (Ai 2,ri_4). Those 



we do not fit, but set them to the values suggested in Ref. [16] using the physical value of ro as 
determined from the fit. This appears to be justified, since we are able to describe our measured 
finite size effects. 

When including NNLO terms into the fit there are four additional parameters to be determined. 
We cannot fit them to the data, because the fits become unstable. In order to be still able to estimate 
systematic uncertainties from NNLO contributions, we include priors for roAi 2, into the fit 

when NNLO terms are included. As priors we use for kM,F = ± 1 and for roAi,2 the values given 



in Ref. [jlj]. 

Our fit procedure can be viewed as first extrapolating the data to the continuum limit and use 
continuum chiral perturbation theory afterwards for the chiral and infinite volume extrapolations. 
For this reason we do not expect any influence of the neutral pseudo scalar meson on the finite 
size effects: in the continuum all three pseudo scalar mesons are degenerate. Note that setting 
D„,p5 = corresponds to a constant continuum extrapolation. Using the boostrap method to 
estimate the statistical uncertainties, we performed the following fits 

1. Fit A: NLO continuum ;tPT, T^j^^ = 0, Anps,/ps = 0. ensembles ^1,2,3,4,6 and Ci,2,3,5 
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2. Fit B: NLO continuum xPT, T^*^ = 0, D,nps,M ^tte^, ensembles 5^2,3,4,6 and Ci ,2,3,5 

3. Fit C: NNLO continuum ^PT, Dmp^j^^ = 0, ensembles ^1,2,3,4,6 and Ci,2,3,5 

4. Fit D: like Fit A, but ensembles B5 and C4 added 

5. Fit E: like Fit A, but ensemble Bj added 

Using the fitted parameters we can then determine low energy constants like £3 4, the chiral con- 
densate £ and the pseudo scalar decay constant in the chiral limit /q. 

Discussion 

The fit results are summarised in table ^. In terms of ;t;^/dof the Fits A, B, C and E provide a 
good description of the data with ;^^/dof ~ 1, whereas Fit D, which includes simulation points with 
mps 600 MeV, has significantly larger ;^^/dof. We conclude from this that ;t^PT is not applicable 
for values of mps > 500 MeV. 

To the contrary, including ensemble Bj as in Fit E, and hence extending the fit-range to a value 
of mps ~ 265 MeV reveals completely consistent results with Fit A. This result makes us confident 
that the extrapolation to the physical point is trustworthy. 

Including lattice artifacts in the fit (Fit B) does change rather little as compared to Fit A, and 
the coefficients D„,pj,/pj are compatible with zero, while the value of X'^/^of is not significantly 
reduced: the differences between the results at j3 = 3.9 and jS = 4.05 can be explained with the 
variance observed in ro/a and Zp. Hence, we have to reduce our (already small) statistical errors 
even further to resolve lattice artifacts in /ps and mps, indicating small lattice artifacts in those 
two quantities. When NNLO terms are included in the Fit (Fit C), the most significant difference 
compared to Fit A is observed for £3. Though this effect is not significant, we include it as a 
systematic error in our final results. 

These findings are visualised in figures l(b' , 2(a) and 2(b) . In figure l(b)| we plot ro/ps as 



a function of the renormalised quark mass ro/^R comparing Fits A and C. In the range 0.04 < 
fol^R < 0.12 the two fits agree remarkably well, while for roHn > 0.12 both fail to describe the 
data. Note that we might be seeing lattice artifacts of the order a^jJ.^ at these large masses, which 
would explain the difference between the results at jS = 3.9 and j8 = 4.05. Similar conclusions can 
be drawn from figure 2(a)| , where we plot (ro/wps)^/ (ro/is) as a function of the renormalised quark 



mass. 



In figure [2(b)| we show the result of Fit B for ro/ps. The three curves correspond to the fitted 
curve at j8 = 3.9 (red), the fitted curve at j3 = 4.05 (blue) and to the continuum curve (black). The 
differences between the three curves are rather small, reflecting the result that ^mpsjps are zero 
within errors. 



3. Conclusion and Outlook 

We have presented an update of the ETMC results for fps and mps and their continuum, ther- 
modynamic and chiral extrapolations. The main difference to the previous analysis in Ref. |^ ^ 
is a new simulation point at j3 =3.9 and the inclusion of ro/a and Zp data into the fit. The 
main results are summarised with £3 = 3.42(8)(I0)(27), £4 = 4.59(4)(2)(13), r'^(2GeV) = 
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Table 3: Summary of fit results. 

(-267 (2) (9) (4) MeV)3 and fj^/fo = 1.0743(7)(3)(4). The first error is statistical, tiie second 
estimates residual lattice artifacts and the third effects from NNLO ;tPT. In addition we have pre- 
sented results indicating that flavour breaking effects are zero within the statistical accuracy, with 
the exception of the neutral pseudo scalar meson mass. 

We thank all members of ETMC for the most enjoyable collaboration. This work has been 
supported in part by the DFG Sonderforschungsbereich/ Transregio SFB/TR9-03. 
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